Angular gaps in radial DLA: two fractal dimensions and non-transient deviations from 

linear self-similarity 
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When suitably rescaled, the distribution of the angular gaps between branches of off-lattice radial 
DLA is shown to approach a size-independent limit. The power- law expected from an asymptotic 
fractal dimension D — 1.71 arises only for very small angular gaps, which occur only for clusters 
significantly larger than M = 10 6 particles. Intermediate size gaps exhibit an effective dimension 
around 1.67, even for M — > oo. They dominate the distribution for clusters with M < 10 6 . The 
largest gap approaches a finite limit extremely slowly, with a correction of order M -0 ' 17 . 



PACS numbers: 61.43.Hv, 05.45. Df, 47.53.+n 

For years, a considerable effort has been devoted to 
the numerical analysis of diffusion limited aggregation 
(DLA), § a model of many natural fractal growth pro- 
cesses. Even though DLA has no characteristic length 
scale, simulations exhibit nontrivial and ambiguous char- 
acteristics. HH In the plane, the fractal dimension D = 
1.71 is usually associated with radial DLA, [1 and D = 
1.66 is associated with cylindrical DLA. J5|-| 

Although the mass of radial DLA obeys the relation 
M ~ Rg' 11 {Rg is the radius of gyration) for all R g 's, 
its geometry varies with R g . Figure |l| compares two 
DLA clusters: One recognizes about AT = 5 main arms 
for M = 10 5 particles, |Sj but TV = 8 or 9 for M = 10 8 . 
H This makes the M = 10 8 cluster look more dense 
(like a filled circular disc). [^]|| However, different in- 
terpretations of TV lead to contradictory results. Lack- 
ing meaningful ways to analyze data, there have been 
very few publications on numerical analysis of DLA since 
1995. Here we present a new systematic description 
that incorporates the arms in a broad context: we study 
the statistics of all the angular gaps between branches. 
While Ref. H pioneered such gap analysis, it had in- 
sufficient data and no data collapse. Application of the 
"e— neighborhood" analysis related the increase in TV 
(and in the apparent density) to an increasing lacunar- 
ity, which could remain consistent with M ~ Rg 71 . The 
two scenarios proposed there - long transient and infinite 
drift - should be replaced by a subtle new one, described 
here. 

Our principal new finding is that, under appropriate 
scaling, and apart from finite size tails, the distributions 
of gap sizes for different M's collapse onto a single func- 
tion, which should also characterize most of the gaps in 
the limit M — » oo. Intermediate size gaps exhibit an ef- 
fective dimension around 1.67, and one needs very small 
gaps, hence very large clusters, to observe the widely an- 
ticipated asymptotic behavior with the unique dimension 
D = 1.71. Interestingly, the two widely quoted values of 
D are both present in radial DLA. 



The intermediate dimension D ~ 1.67 is not a small 
sample correction to asymptotic self-similarity. Indeed, 
the number of gaps in this zone changes little with M, 
and they remain observable for all M . Therefore, DLA 
cannot be modeled by a unique fractal dimension, that 
would correspond to a fractal that is invariant under lin- 
ear transformations. For many phenomena, like perco- 
lation, random forms of such "linear fractals" are good 
models and non-random linear fractals provide illuminat- 
ing "cartoons". In contrast, DLA (and probably other 
phenomena) require steps beyond linear self-similarity. 
Our procedure suggests such a quantitative step. 

The data collapse of the gap distribution implies that 
typically, the largest angular gap (which gives some mea- 
sure of TV) approaches a finite limit, which seems to be 
non-zero. However, there is an amazingly slow correc- 
tion, of order M"T, where 7 = (2 - D)/D = 0.17. This 
correction may or may not be related to the corrections 
recently found using a conformal map approach. pT| , ^2| 

We performed large scale simulations of off-lattice ra- 
dial DLA, and focused on the gaps on circular cross-cuts. 
An earlier analysis, |jj based on M < 10 7 , yielded an ap- 
parent fractal dimension D = 1.65 ± 0.01, which is not 
consistent with the widely accepted D — 1.71. Our soft- 
ware is as described in Ref. |l3) , but uses a single diffus- 
ing particle each time, as in the definition of DLA. jij We 
generated N C (M) — 1000 clusters for each of the 'masses' 
M = 10 5 , 10 5 - 5 , 10 6 , 10 6 - 5 and 10 7 , plus N C {M) = 60 
additional clusters for M = 10 7 5 and for 10 8 . 

After growing several DLA clusters, we endeavored to 
keep to the edge of the non-growing region. Trial and 
error ||J| singled out particles whose distance R from 
the origin is within 0.75i? g - 1 < R < 0.75R g + 1 (in 
units of the particle diameter). The gap 9 is the dif- 
ference in the polar angles between consecutive particle 
centers. To avoid "gaps" between touching particles, we 
keep only the values > l/0.75i? 9 , larger than one parti- 
cle diameter. We then pool the gaps from all the available 
clusters of size M into a single list, which is sorted from 
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largest to smallest. Thus each gap is characterized by its 
size and by an integer r that represents its rank in the 
sorted list. Here we use the mean-rank, p = r/N c (M), 
which becomes an A/— independent continuous variable 
for N c — > oo. 

For an ideal linearly self-similar fractal of dimension 
D, one expects a power law, 9 oc p~ a , where a = l/D c 
and D c = D — 1 is the fractal dimension of the intersec- 
tion. Q In contrast, Fig. g shows that 6(p,M = 10 s ), 
based on the available iV c (10 8 ) = 60 clusters, exhibits 
four distinct regions: Region I, with p < 10, reflects the 
statistics of the largest gaps. As in the Zipf- Mandelbrot 
distribution, |Q the flattening of the curve in this re- 
gion is a token of an upper cutoff (though not necessar- 
ily a sharp one) on the size of gaps. Both Regions II 
and III appear to be linear (see dashed lines in Fig. ||). 
In Region II, with 10 < p < 50, the apparent slope is 
an = 1.5, corresponding to Du = 1.67. In Region III, 
with 50 < p < p m a X /3, where p max (M = 10 8 ) = 1890 
is the maximal p for which 9 is defined (i.e. the average 
total number of gaps), the apparent slope is am = 1.33, 
corresponding to Dm = 1.75. Region IV, with p ma x/3 < 
P < Pmax, reflects finite-size corrections: as M — > oo one 
has pmax — * oo, and (as far as we can tell) Region III 
extends to infinity, reflecting the asymptotic behavior of 
DLA. Since p ma x increases with M, it also follows that for 
M smaller than about 10 6 , Region III merges with Re- 
gion IV, and one would mistakenly identify D from the 
higher apparent slope of Region II. This clearly demon- 
strates why one must have very large clusters to make 
statements about the asymptotic behavior of DLA. 

Visually, Region I defines the arms. As M — > oo, Re- 
gion III collapses into the observed thickening 'trunks'. 
Region II determines the structure of the clearly separate 
intermediate branches, which exist for all M. 

A quantitative way to describe the crossover between 
Region II and III would involve corrections to the asymp- 
totic power law. Lacking enough data, we represent 
the data we have by an ad-hoc parametric power ex- 
pansion, 6{p) = J2k=T a kP~ ka = Ylk=i a kX k , where 
x = p~ a ■ Using a = 1.4, this form ensures the asymp- 
totic dimension D = 1.71 for large p. A least squares 
fit with A: max = 6 to the data for 10 < p < p max /3 at 
M = 10 s yields a k = 3.702, -126.7, 1.566 x 10 4 , -8.045 x 
10 5 , 1.901 x 10 7 , -1.694 x 10 8 . The approximate effec- 
tive dimension, -D c ff(p) = 1 + l/a c ff, where a e g(p) = 
-<2(log#)/(f(logp) = af J|, is shown in Fig. [| Indeed, 
D e s reaches a minimum of about 1.66 close to p = 20 
(in Region II), and a maximum of about 1.74 around 
p = 80 (in the beginning of Region III). It approaches 
the asymptotic D — 1.71 at larger p. 

To demonstrate scaling and data collapse, we scale 9 

as 



S(M) = / 8(p,M)dp. 



6(p,M) = S(M)e n (p), 



(1) 



(2) 



S(M) represents the largest gaps. G should be large 
enough in order to filter more of the statistical variance, 
but not too large as to include finite-size corrections. The 
results are not very sensitive to G (see below) . The choice 
G = 50, at the beginning of the asymptotic slope (Re- 
gion III), gives good results and facilitates some of the 
theoretical analysis presented later. Figure [I| shows that 
9 n (p,M) is indeed independent of M, except for finite- 
size correction deviations. 

The curves in Fig. f| end at p max (Af). The figure 
clearly shows fixed distances between consecutive end- 
points, indicating power-law dependence of both p max 
and of #(pmax) on M. Since p max is the average number 
of gaps (or particles) on the annulus of radius _R, we 



expect p max oc R 



D-l 



Mi - 1 )/ . The smallest gap 



is then expected to scale as l/R oc l/M 1 / . We thus 
expect the small gaps to scale as 



8(p, M) = M- 1 ' fip^^-V/M 1 ' ). 



(3) 



Indeed, Fig. || shows M X I D Q plotted vs. p 
for various M's. There is a clear data collapse on the 
right-hand side, in agreement with Eq. (^). 

We next use the collapse of the normalized gaps in Fig. 
|] to infer the asymptotic distribution of the original (non- 
normalized) gaps as a function of M. As expected for p < 
Pmax/3, n depends only on p, and not on M. The in- 
dependence of the original gaps, 9(p,M) = S(M)9 n (p), 
is thus determined by the scale factor S(M). Integrating 
over p on both sides of Eq. ([j]), the left hand side be- 
comes J Pmax 0(p, M) = 2tt — X, equal to the average sum 
(over all cluster realizations) of all the gaps; X is the sum 
over the gaps that are smaller than one particle diameter. 
We expect X oc p max /R g , and therefore X = i?iAf~ 7 , 
with 7 = (2 - D)/D = 0.17. On the right hand side, 
we divide the integral into three parts: < p < G, 
G < p < p m ax/3 and p max /3 < p < p ma x- Equation 
(||) implies that J Q 9 n dp = 1. The second integral be- 
longs to Region III, where we write 9 n (p) ~ Aip~ a , and 
Ax = ai/S^lO 8 ) = 0.77 ± 0.01 is independent of M (as 
seen from Fig. |). This yields B 2 \G~ (i - (p max /3)-% 
with (3 = (2 - D)/(D - 1) and B 2 = Ai/(3. Since 

Pmax OC 

M (D-1)/D 

we conclude that the integral over 
the right hand side is equal to S(M)(A — S3 A/ -7 ), where 
A = I + B2G- 13 = 1.40±0.01 and B 3 is a constant (which 
is also affected by Jp max / 3 9 n dp ~ Af ~ 7 , due to Eq. (||)). 
Our conclusion is that for large M one has 



S(M) 



2tt - BiAf- 7 
A - B3M- 7 '■ 



(4) 



and hence 6»(p, M = 00) = (2n/A)9 n (p) = 4.50 n (p). 

While the distribution of gaps does converge to a limit 
for Af — > 00, this convergence is amazingly slow, being 
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governed by the small exponent 7: The correction term 
in Eq. (0) reduces only by a factor of I0 3x017 ~ 3.2, 
between M = 10 5 and M — 10 s ! This means that great 
caution should be taken when extrapolating various mea- 
surements in DLA, and that the asymptotic limit is at- 
tained only for extremely large M . 

We next examine the statistics of the largest gaps, 
#max- For a linearly self-similar fractal, one expects to 
find the same distribution of the # max 's for different clus- 
ter sizes. # m ax ~ 2ir/Af may indicate that the number 
of main branches is around J\f. It has been shown that 
#max decays slowly with M, |3| which means that TV in- 
creases, raising doubts on the self-similarity of DLA at 
intermediate values of M. 

If, instead of pooling, we had studied the gaps of a 
single cluster, then we would have # max = 9(r = 1). 
However, the mean-rank p = 1 has a different meaning: 
9(p = 1, M) is the Y c 'th largest gap out of all the gaps 
taken from N c clusters. Since the graph of 8 n (p,M) is 
relatively flat near p = 1, we still expect 6(p = 1, M) to 
represent the order of magnitude of the average largest 
gap. Reading 9 n (p = 1) = 0.1 from the universal curve in 
Fig. |, we estimate 9(p = l,M = 00) = (2ir/A)6 n (l) ^ 
0.45, in radians, yielding an asymptotic average maximal 
gap around 26° and implying a finite number of arms. 

The histogram of the 6> max 's for different clusters with 
M = 10 7 is shown in Fig. || It is quite broad. A fit to 
a log-normal form (also shown) yields p = (In # max ) — 
-0.50 and a = ((m0 max - pf) 0.22. We find a very 
slow decrease of p with M, approaching a finite limiting 
value of the order of magnitude given above. 

The broad distribution of # max implies that under pool- 
ing, the largest pooled gaps spread over a broad range 
of p. Pooling convolutes the distributions of two non- 
independent quantities: the gaps for each cluster and the 
different clusters. For additional insight into these distri- 
butions, we also considered the correlations between gaps 
within each cluster. We find several interesting results: 
the average correlation between the largest gap and the 
r'th gap for M — I0 7 turns out to decrease from 1 to 
(around r = 5), and to approach a constant around —0.5 
for r > 10. Negative correlations were expected, since 
the total sum of the gaps must equal 2ir. 

In collaboration of Eugene Vilensky of Yale, we also 
tested our results' sensitivity to the special choice of R ~ 
0.75i? g . An increase of R beyond this value first leads to 
little change (confirming then to opposite changes 

in different regions: Very large gaps come to dominate in 
Region I, moving the corresponding "flat" portion of the 
curve to higher values of 9. The remaining gaps (other 
than the few largest) add to a decreasing total angle, but 
with little change in distribution. 

In conclusion, we have identified robust structural fea- 
tures of radial DLA, which act separately in character- 
izing the arms, the intermediate and the fine branches. 
Although all fall within a single scaling function, the non- 



trivial gap distribution function reflects deviations from 
simple linear fractality. Since these deviations character- 
ize the radial geometry, it would be interesting to check 
if they also appear in the (possibly simpler) cylindrical 
case. This is particularly intriguing, since both involve 
D w 1.66. Finally, since our data 'only' extend up to 
M = 10 8 , there is no guarantee that this is the end of the 
story. [^5) Future work should tell if larger scales bring 
more surprise. 

Henry Kaufman kindly adapted his DLA program to 
serve our needs. This work has been supported in part 
by the German-Israeli Foundation (GIF). 
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FIG. 1. Two radial DLA clusters, with M = 10 5 (top) and 
M = 10 8 (bottom) particles. 
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FIG. 2. The gap size 6 vs. the mean-rank p. The graph is 
based on 60 clusters with 10 s particles. The dotted vertical 
lines separate the four regions discussed in the text. Dashed 
straight lines run through Regions II and III. 
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FIG. 3. Effective dimension _D of r(p). 
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FIG. 4. Data collapse of 8 n (p,M). The numbers indicate 
the cluster size M. 
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FIG. 5. M 1/D 6 vs. p 1/(I5 - 1) /M 1/I5 for various values of M. 
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FIG. 6. Histogram of the maximal gap (in radians) in each 
cluster, for the 1000 available clusters with M = 10 7 . 



